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Abstract 

An optimization method used in image-processing (metamorphosis) is found to imply Euler's equa- 
tions for incompressible flow of an inviscid fluid, without requiring that the Lagrangian particle 
labels exactly follow the flow lines of the Eulerian velocity vector field. Thus, an optimal control 
problem and an optimization problem for incompressible ideal fluid flow both yield the same Euler 
fluid equations, although their Lagrangian parcel dynamics are different. This is a result of the 
gauge freedom in the definition of the fluid pressure for an incompressible flow, in combination with 
the symmetry of fluid dynamics under relabeling of their Lagrangian coordinates. Similar ideas 
are also illustrated for SO(N) rigid body motion. 
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Hamilton's principle for ideal fluid flow might be summarized by saying that water moves 
as well as possible to get out of its own way jl]. This phrase could make the challenge of 
optimal control of fluids seem daunting, particularly when combined with Le Chatelier's 
principle that any complex system will respond to resist the effects of an external force. The 
question pursued here is whether Euler's fluid equations (EFE) represent optimal control, 
or only optimization. As it turns out, the geodesic flow represented by the EFE is found to 
arise from either formulation. 

Definition 1 (Optimal control) 

Optimal control problems consist of Q]: 

• A differentiable manifold M on which state variables Q G M evolve in time t during 
an interval / = [0, T] along a curve Q : I — ► M from Q(0) = Qo to Q(T) = Qt, with 
specified values Qo, Qt G M; 

• A vector space V of control variables U G V whose time dependence U : I — * V is at 
our disposal to affect the evolution Q(t) of the state variables; 

• A system of first-order evolutionary partial differential equations (state equations) 
defined on the tangent bundle TM by the vector field F:Mx^-> TM, 

Q = F(Q,U), (1) 

and introduced as a constraint that relates the unknown state and control variables 
(Q(t),U(t)) : I -f M x V, and 

• A cost functional depending on the state and control variables 

S:= [ T £(Q,U)dt, (2) 
Jo 

whose minimization is the goal, subject to the prescribed initial and final conditions at 
Q(0) and Q(T), and the state equations ([1]). The integrand L : M x V — > R is called 
the Lagrangian, which is assumed to be continuous and continuously differentiable on 
M x V. 



Thus, solving a standard optimal control problem requires finding time-dependent state and 
control variables (Q(t),U(t)) that optimize a given cost functional S subject to the exact 
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enforcement of an evolutionary system of state equations (JTJ), while satisfying prescribed 
initial and final conditions, Qo and Qt- 

The optimal control problem is traditionally formulated by introducing a pairing ( ■ , ■ ) : 



gent space of costate (Q,P) G T*M variables. In terms of this pairing, the state equation 
may be enforced as a constraint on variations of the cost function, by using the classical 
method of Lagrange multipliers. Thus, the optimal control problem in this formulation 
becomes, 



for which standard methods of variational calculus are available. (The variations are to be 
taken at fixed values of the endpoints Qo and Qt-) 

Definition 2 (Optimization by metamorphosis) 

The term metamorphosis refers to a class of optimization methods used for performing image 
registration by finding the optimal flow along a curve in the group of of diffeomorphisms 
Diff(M) (smooth invertible maps with smooth inverses) acting on a differentiable manifold 
M of image properties (states) Q G M defined over a given spatial domain T> Q . Such an 
optimal flow is sought as a geodesic time-dependent curve with respect to a certain metric 
on the tangent space TDiff(M). Hence, one chooses control variables U G V — TDiff(M). 
One version of optimization by metamorphosis replaces the optimal control problem by 
the minimization of a sum of norms integrated over the control time interval. This version 
of metamorphosis may be expressed as the following optimization problem, 



for a positive real constant a 2 and a chosen norm || • || : TDiff(M) — > R, while satisfying 
prescribed initial and final conditions, Qo and Qt- 

The optimization problem (j3J) is solved by finding controls U (t) that steer the state variable 
Q along a time-dependent curve Q(t) leading from Qo to Qt, obtained by minimizing the 
cost of applying the controls S in (J2J), while enforcing the state equations ([1]) within a certain 
tolerance a. The solutions (Q(t), U(t)) of an optimal control problem and its corresponding 
problem of optimization by metamorphosis need not coincide, even when their cost functions 
and state equations are the same. 



TM* x TM — > R between the tangent space of state variables (Q, Q) G TM and the cotan 




(3) 
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Optimization by metamorphosis finds applications in problems of image registration and 
recognition of patterns, e.g., in the analysis of medical images obtained using MRI, CT 
and other imaging technologies [3, 0|. In these applications, shapes must be matched, or 
at least compared to one another. The matching or comparison procedure is formulated 
as an optimization problem whose goal is to minimize the sum of the chosen norm on 
the tangent space of the path plus another norm associated with the estimated error, or 
tolerance, of the measurement process. This formulation of optimization by metamorphosis 
has the advantage of introducing a Riemannian structure that ensures that the extremals in 
the distance between the images are genuine minima, rather than being saddles. For more 
discussion of the metamorphosis approach, see [3|, |4j. 

In addition to image analysis, an interest in controlling diffeomorphisms also arises in ideal 
fluid flows. When applied to fluid dynamics, optimization by metamorphosis introduces a 
penalty defined by a metric on the tangent space of the inverse flow (also known as the 
back-to-labels map for fluids). The penalty introduces an additional cost in the kinetic 
energy of labels whose paths deviate rapidly from the Eulerian characteristics following the 
forward flow. However, the metamorphosis approach does not constrain the fluid labels to 
follow exactly along the forward flow lines. 

Objective The present paper shows that optimization by metamorphosis recovers the clas- 
sical Euler fluid equations (EFE), even though the metamorphosis approach does not require 
that the fluid labels exactly follow the characteristic curves of the Eulerian fluid velocity. 
This result implies that the standard Lagrangian representation of fluid dynamics as la- 
belled fluid "parcels" that are carried along characteristic curves of the Eulerian velocity is 
sufficient for deriving the EFE; but it is not necessary. 

EFE for incompressible inviscid flow were identified as an optimal control system by deriving 



them using the Hamilton-Pontryagin principle in [5] . This derivation of EFE from an optimal 
control problem recovered Arnold's interpretation of EFE as geodesic flow on the volume- 
preserving diffeomorphisms It also recovered a standard representation of the Euler 
equations called the impulse equations [7, 8, 9]. Applications of the impulse representation 
in numerical simulations of EFE were discussed in [loj]. One disadvantage of the optimal- 



control formulation of EFE in 



5j was that it enforced the pointwise physical constraint 



that Lagrangian fluid parcels are frozen into the flow. That is, the Lagrangian particles 
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were required to follow exactly along the flow lines of the Eulerian velocity vector field. This 
strict pointwise constraint led to a variational principle whose extremals were not necessarily 
true minima. Instead, they could have been saddles; so the controllability of Euler fluid flows 
was left as an open question. Controllability of Euler fluid flows will also be a moot point 
in the present work. Geometric control |J of Euler fluids using Lie methods is reviewed 
expertly in 11]. 

In the present paper, we first derive EFE from an optimal control problem using the Clebsch 
variational approach that enforces the frozen-in particle constraint via the back-to-labels, or 
inverse map, as done previously in for the forward map. This approach produces a sym- 
metric form of EFE analogous to the symmetric, double-bracket form of the iV-dimensional 
rigid body dynamics found in [l2]. As a result, the Clebsch representation of the EFE 
may be written as a coupled system of double-bracket equations in analogy to the corre- 
sponding representation of the SO(N) rigid body. We then re-derive EFE using the method 
of metamorphosis. From the viewpoint of fluid dynamics, it is interesting that the same 
fundamental EFE appear, even though the metamorphosis approach does not impose the 
strict requirement that Lagrangian labels follow along characteristic curves of the Eulerian 
velocity vector field. 

Plan The paper has four parts. The first part reviews the properties of the optimal control 



problem for the SO(N) rigid body 12|, 



13 



14j . These properties include a Q <-» P exchange- 



symmetric canonical Hamiltonian formulation that may be rewritten as a coupled system 
of double-bracket equations on SO(N) x SO(N). The remarkable features of the SO(N) 
rigid body provide a model for the paper's subsequent development. The second part shows 
that an optimization problem for the SO(N) rigid body leads to a derivation of the same 
dynamical equations as found from its optimal control problem. The third part derives 
the EFE from an optimal control problem that constrains Lagrangian trajectories to follow 
Eulerian flow characteristics locally at every point. This approach results in the well-known 
Clebsch formulation of the EFE [jjj]. In the Clebsch formulation, the same properties 
of canonical exchange symmetry and double-bracket dynamics emerge as for the SO(N) 
rigid body. The fourth part derives the EFE from an optimal control problem based on 
metamorphosis that constrains the Lagrangian trajectories only in an L 2 sense. The EFE 
still reappear from this optimal control problem, albeit with a modified pressure and broken 
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exchange symmetry in the canonical equations. The two differences are immaterial, though, 
because the term that breaks exchange symmetry in the metamorphosis approach appears 
in its resulting EFE as merely a redefinition of pressure, i.e., a gauge transformation, which 
has no effect at all on the solutions of the Euler fluid equations. 

Optimal control of the iV-dimensional rigid body Let us recall the standard and 
symmetric forms of the equations for a rigid body in N dimensions. We first recall that the 



left- invariant generalized rigid body equations on SO(iV) may be written as 16 

Q = Qn, M = [M, fi] , (5) 

where Q G SO(iV) denotes the configuration space variable (the orientation of the body), 
Q = Q~ X Q G so(N) is the body angular velocity and 

M := J(Q) = Af2 + f2A G so(N)* (6) 

is the body angular momentum. Here the positive definite operator J : so(N) — > so(N)* is 
symmetric with respect to the matrix trace inner product 

(A,B)= 1 -tr(A T B). (7) 

The diagonal matrix A satisfies A, + Aj > for all i ^ j. For n = 3 the elements of A$ are 
related to the standard diagonal moment of inertia tensor I by l\ = A 2 + A 3 , J 2 = A 3 + A 1; 
J 3 = Ai + A 2 . 

The equation M = [M,Q] is readily checked to be the Euler-Poincare equation on 5o(N)* 
for the Lagrangian 

/(fi) = i(fi,J(fi)). 

Theorem 1 The left-invariant rigid body dynamics is given by the symmetric system of 
first-order equations 

q = Qn, P = pn, (8) 

where Q is regarded as a function of Q and P via the equations 

n := J~\M) eso(N), M := Q T P — P T Q. (9) 

Sketch of proof: The result follows by direct substitution of equations © into (jSJ), as 
shown in 3, Q, 3- ■ 
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Remark 1 The exchange-symmetric system (jHJ) for SO (N) rigid body motion is canonically 
Hamiltonian for H(Q, P) = \{J~ l M(Q, P), M(Q, P)) with Poisson bracket {Q, P} = Id. 
The map M : T*SO(N) — > so(N)* in ([9]) is the cotangent-lift momentum map for the (left) 
action of the Lie group SO(N) on itself. Substituting equations ( 
them on SO(N) x SO(N) as a coupled double-bracket system \l£, 



9]) into (jHJ) reformulates 

3,Qy. 



The symmetric form of the n-dimensional rigid body equations on SO(N) in (jHJ) is interesting 
for our purposes here, because these rigid-body equations can also be derived from the 
following optimal control problem. 

Definition 3 (Optimal control problem for the SO(iV) rigid body) 

Let T > 0, and let Q , Qt G SO(iV) be fixed. The rigid-body optimal control problem is given 

by 

T 

(10) 



(11) 



min - (U,J(U))dt, 

UGso(N) 4 J 


subject to the constraint on U that there be a curve Q(t) G SO(iV) such that 

Q = QU Q(0) = Q , Q(T)=Q T . 



In the framework of the symmetric representation of the rigid body equations (jHJ) the fol- 
lowing theorem may be proved. 



Theorem 2 (Bloch et al. jl2j |) The rigid body optimal control problem given in Defini- 
tion [3] has extremal evolution equations (jHJ) where P is the costate vector given by the Pon- 
tryagin maximum principle. The optimal control relation in this case is given by 



U = r\Q T P-P T Q) 



(12) 



Definition 4 (Optimization problem for the SO(iV) rigid body) 

Let T > 0, and let Qo,Qt G SO(iV) be fixed. The rigid-body optimization problem is given 

by 



1 

mm - 

U£so{N) 2 



{U, J(U)) + -\\Q - QU\\< 



0~~ 



dt, 



for a positive real constant a 2 , a given metric || • || : TSO(iV) 
endpoint conditions on Q that the curve Q(t) G SO(iV) satisfy 



(13) 

and subject to the 



Q(0) = Qo, Q{T) = Q T 



(14) 
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Remark 2 The problem statement (Q2|) optimizes the cost / fiOj) 0/ «sm<? f/ie controls, for a 
given tolerance (a) in satisfying the state equations. 

Suppose the metric || ■ || : TSO(iV) — > IR in (TT3T) is given by 

\\ Q -QU\\ 2 = ti {{Q -QU) T K{Q -QU)^j , with K T = K , (15) 

as obtained from the trace inner product for matrices in (171), where, as expected for a metric, 
K is symmetric. In this framework, the following theorem may be proved. 

Theorem 3 (Optimization equations for the SO(iV) rigid body) The rigid body op- 
timization problem given in Definition @] yields extremal evolution equations, 

K(Q - QU) = o 2 P , p = -PU T = PU, (16) 

obtained from variations of the cost function ( TTBl in the variables Q and Q, respectively. 
The optimal controls in this case are again given by equation CQ| ). 



Sketch of proof: The extremal evolution equations in ( |T6l) follow from stationarity of 
the cost function ( |T3l) under variations in Q and Q, respectively. Independently requiring 
stationarity of the cost (fTBl) under variations of the control variable U yields 

J{U) = \{Q T P-P T Q). (17) 

This reproduces the momentum map in (Q and thereby recovers the optimal control relation 
in (frj). ■ 

Theorem 4 (SO(iV) rigid body motion optimizes ( 1131) ) For any value of the toler- 
ance a 2 , the optimization problem for the SO(iV) rigid body in Definition® yields an evolu- 
tion equation in the same form as the angular momentum equations for the rigid body, when 
written as 

M = [M, U] , with M := J[U) , (18) 



and J(U) is given by the control relation [T7\ ). 



Sketch of proof: The extremal evolution equations ( 1161) combine with the optimal control 
relation ( 1121) to produce 

d -^l = [J(U),U] + yi^ 1 - K- T )P , (19) 



Deriving Euler's equations 



D. D. Holm 



9 



The last term vanishes, by symmetry of K, thereby recovering the motion equation (fl8|) for 
any value of a 2 . ■ 

Next, we will derive the symmetric form of the Euler fluid equations on the volume-preserving 
diffeomorphism group SDiff. These equations correspond to the symmetric equations (JE]) on 
SO(N). We will also compare the equations that result for the control dynamics as they 
are obtained from the optimal control, and optimization approaches. 

Optimal control for fluids: Clebsch approach 

Let us begin by explaining the statement made in the Introduction that metamorphosis 
penalizes fluid parcels that deviate from their flow lines, but does not force them to follow 
exactly along the flow lines of the fluid velocity. We recall that Lagrangian fluid dynamics 
provides evolution equations for particles moving with a fluid flow. This is typically expressed 
via the forward flow map g t := g(t) from the fluid reference configuration at time t = (for 
which go = Id) to the fluid flow domain T> at a later time t. The flow lines are given by 
Xt = gtl with xq = I, the vector label. (Here vector and covector indices are understood, 
but not written explicitly.) The forward map is taken to be a diffeomorphism, so the 
flow preserves topology. That is, the fluid particles are imagined to be unable to cavitate, 
superimpose or jump. 

For an n- dimensional fluid flow, the flow line relation x t = gtl of the flow map g t : W 1 x K i— ► 
R n specifies the spatial position at time t of the fluid particle that has label I = xq. The 
inverse map g^ 1 is the back-to-labels map, {23}], which gives the label of the particle that 
occupies position x at time t as the vector function l t = g^~ x. The Eulerian velocity field 
u t = u(x,t) gives the velocity of the fluid particle that occupies position x at time t as 

9t = u t og t or x(l,t) = u(x(l,t),t). 

The vector components of particles whose labels l t are frozen into an ideal fluid flow each 
satisfy the advection law obtained from the time derivative (tangent) of the back-to-labels 
map l t = g^ x = l{t,x), 

— : =d t l + u-Vl = 0. (20) 



The quantity v = —dtl(t,x) = u ■ VI is called the convective velocity 18|. Satisfying 



v = Adg-iu, the convective fluid velocity v is to the Eulerian spatial fluid velocity u as body 



angular velocity is to spatial angular velocity for a rigid body 



24) 
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The optimal control problem for EFE may be posed using the vector L 2 pairing (a , b) 
§ v a T b d n x as 



T 



min / £(w)dt with £(u) = -(u,u), (21) 
«*(•) J 2 
o 

subject to divu = 0, d t g — u o g, with g(l, 0) = go(l), g{l, T) = gr(l) fixed and, for flow in 
all of space, suitable conditions at infinity. 

The optimal control problem (]2ip is of course identical to the standard Hamilton principle 
for ideal fluid mechanics and it has been solved previously in terms of the forward map 5 ] . 
However, that approach using the forward map did not produce the symmetric form of the 
EFE. Here, we solve problem (12T1) and find the symmetric form of the EFE by introducing 
Lagrange multipliers (n, k) into the kinetic-energy cost that impose incompressibility and 
the pointwise tangent relation (TSUI) along the inverse flow, i.e., the back-to-labels map l(t, x), 



T 

S(u, I, n, k) = J (i(u) + (n, d t l + u- VI) - jfc, div u)^ )dt. (22) 
Constraint Constraint 

The problem may be recast as: minS", subject to the two constraints divw = and dtl + 
u ■ VI = 0, as well as spatial boundary conditions and endpoint conditions that Z(0,x) and 
1{T, x) are fixed. 

Remark 3 The absence of weak solutions in L 2 for the Euler equations is a well-known 
problem, both for the initial value problem and for the geodesic minimization problem, see, 
e.g. 



19i . However, for the stronger norms used in typical applications of metamorphosis 
in image analysis Jj/ and with the H 1 norm for the averaged Euler equations \2w , one may 
assume that all entities exist, functions are smooth and integrals are finite. Consequently, 
we will perform the calculations below with a sufficiently smooth Lagrangian £(u) that these 
assumptions hold in either the optimal control problem in equation / f^j) or the optimization 
problem in equation (2l\ ). We will then specialize the results to the case l[21\) ofu&L 2 when 
referring to the Euler equations, knowing that the result is only formal in this case. 



For the optimal control problem in equation (|22|) . the following results may be proved for a 
sufficiently smooth Lagrangian £(u). 
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Theorem 5 The extremals of S(u,l,ir,k) are given by 

- + 7 rV/ + gradfc = 0, ^ = = - 



(23) 



and divu = with 



D 

TTt 



d t + u-W. 



Sketch of proof: The variations of S(u, g, n, k) yield 

T 

5S = 

o 



7 — h 7rV/ + grad k,5u) — (5k, div u) 
ou 



+ (Sir, d t l + u ■ VI) - (d t n + div nu, 51) dt + 



ir 



SI 



T 



System (123|) follows immediately, upon using divw = and noting that 8u{po, t) = 51(0, x) = 
5l(T,x) = 0. m 

Remark 4 The (ir,l) equations (123|) are symmetric under exchanging tt <-> I. These are the 
analogs for SDiff of the P Q symmetric equations (jSJ) for SO(N). The representation 
of the control variable u G X (a divergenceless vector field) in terms of the state variables 
7i, I, k is given by inverting the map m(u) G X* in 

5t 

m(u) := — = -TiVl-Vk, (24) 
5u 

obtained from the variations in u of the cost function S. The Clebsch map m : T*SDiff — > X* 
in fT2~4"l) is the cotangent-lift momentum map for the (right) action of SDiff on itself. 

Theorem 6 Substitution of the Lagrangian £(u) in (I2ip into system (|23|) for the extremals 
of S yields EFE. 



Sketch of proof: System 
for ideal fluids 



implies the following version of Kelvin's circulation theorem 



d f S£ . 
— (b — ■ ax 



dt J c{u) 5u 



c(u) L 



D 5£ ,„ , T 51 
mTu + {Vu) 'Tu 



■ dx 



(Dt:.. Ml Dk\ n 

dl + ird h d = 

c(u)\Dt Dt Dt) 



(25) 



whose RHS vanishes upon using (T23]) and noting that the integral of an exact form vanishes 

when taken around the closed loop c(u) moving with the Eulerian velocity u. For the Euler 

case ( 1211) . in which SijSu = u, this calculation recovers EFE as 

Du ( Dk u 2 " 
= -V + — 



Dt 



\Dt 



Vp and div u — . 
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Vorticity dynamics By taking the curl, the Euler equations for the vorticity w = curlw 
are found to be 

dw 

— = [w,u\, w = V7r x V/, divw = 0, (26) 

where [w, u] — w • Vu — u ■ Vw is the Lie bracket of the divergenceless vector fields X on the 
flow domain T>. The interpretation of (1261) is that the symplectic 2-form dir A dl is frozen 
into the forward map. 

Remark 5 When 5£/5u = u the system ( 123]) for the extremals of S recovers the classical 



Clebsch representation of ideal incompressible fluid flow flolj . The vorticity equation ( 1261) is 
canonically Hamiltonian for 

H(7r 1 l)= 1 -((Ar 1 w(7rJ),w(7r 1 l)}, 

with Poisson bracket {tt(x),1(x')} = —5(x — x'). The map w : T*SDiff — * X* on T> is the 
cotangent-lift momentum map for the action of SDiff on itself. Substituting equations ([9]) 



into (E]) reformulates them on SDiff x SDiff as a coupled double-bracket system }2n] 



Metamorphosis and the optimization problem for fluids 

Now we come to the point of formulating the optimization problem ( 12~TT) for EFE whose 
extremals are sure to be minima. This is the metamorphosis formulation, which solves 
problem (121]) by introducing an additional norm into the kinetic-energy cost that imposes 
incompressibility as a constraint and treats the tangent relation (120]) for the back-to-labels 
map as merely an optimization penalty. This is 

T 




Penalty Constraint 

The problem may be recast as: minS 1 , subject to divu = 0, spatial boundary conditions, 
endpoint conditions that 1(0, x) and l(T, x) are fixed and penalize for the error in the L 2 
norm \\dtl + u ■ Vl\\ 2 L2 - For a 2 > 0, when extremals exist they will be minima. 
For this problem, the following results may be proved. 

Theorem 7 The extremals of S(u,l,k) are given by 

5£ Dn Dl „ 

- + ,V i + gr a" = 0, Di = °' Di = **< (28 » 
and divn = with — :— dt + u ■ V. 
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Sketch of proof: After defining the canonical momentum by a 2 ir = d t l + u ■ V/ obtained 
from varying d t l, the other variations of S(u, I, k) yield 

T 

5S = J ((^+irVl + gTedk,6u\-(5k, divu) 
o 

i T 



(d t 7r + div 7ru, 51) dt + 



71 51 







System (1281) follows immediately, upon using divu = and noting that 5u(oo, t) = 51(0, x) = 
5l(T,x)=0. m 

Remark 6 Equations f l28|) break the exchange symmetry tt I, which is only restored in 
the limit a 2 — > 0. 

Theorem 8 Substitution of the Lagrangian t{u) in fl2Tl) into system fl28l) for the minima 
of S again yields EFE, for any value of a 2 > 0. 

Sketch of proof: System ( 1281) implies the following version of Kelvin's circulation theorem 



d f 5£ . 
— (b — ■ ax 



dt J c (u) 3 U Jc{u) 



Dt 5u 5u 



■ dx 



fDiTj, Dl Dk\ n 

d/ + ?rd hd =0, 

c(«)V^ Dt Dt) 

whose RHS vanishes upon using ( 1281) and noting that ixd(Dl / Dt) = d(<T 2 7r 2 /2) is exact, so 
it does not contribute to the integral taken around the closed loop c(u) moving with the 
Eulerian velocity u. For the case 5lj5u = u, this calculation recovers EFE as 

Du ^(Dk u 2 9 7r 2 \ „ 

= - V 1 h a 2 — =: - Vp, 

Dt \Dt 2 2 J 1 ' 

and divw = yields pressure p independently of a 2 . ■ 

Remark 7 One concludes that an optimal control problem and an optimization problem 
for incompressible ideal fluid flow both yield the same geodesic Euler fluid equations for 
their control relations, although they represent different Lagrangian dynamics. The EFE 
equations derived from either the optimal control problem based on the Clebsch approach, 
or the metamorphosis optimization approach are equivalent up to a gauge tranformation of 
the pressure, which is immaterial for incompressible flow, since the pressure in this case 
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is determined independently by preservation of the constraint divu = 0. The emergence 
of equivalent EFE from either optimal control or optimization means that the forward map 
in the Lagrangian picture of ideal fluid dynamics implies the Eulerian picture, but not vice 
versa. 

The label paths defined by the inverse map It = g^ l x may 'forget their way" and their 
dynamics may diverge from the requirements of the Eulerian velocity characteristics as in 
for any finite value ofo~ 2 >0 without changing the optimal outcome. Taking the limit of 
sigma to zero yields exactly the same relations for the label dynamics of both the optimization 
problem and the optimal control problem. However, this result is immaterial to the emergence 
of the EFE, because both problems yield the same Euler fluid equations, even though their 
Lagrangian parcel dynamics are different. The result arises from the gauge freedom in the 
definition of the fluid pressure for an incompressible flow. One notes that the gauge freedom 
in defining the pressure for incompressible fluid flow is independent of the symmetry of fluid 
dynamics under relabeling of their Lagrangian coordinates. The relabeling symmetry implies 
the conservation of circulation in and holds for both Lagrangian cost functions ft2E) and 
(2l\ ), independently of the definition of pressure. If the relabeling symmetry were broken to 
a subgroup corresponding to invariance of fluid properties appearing in the thermodynamic 
definition of pressure, the two formulations may produce a nontrivial difference. This feature 
will be investigated elsewhere. 
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